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Abstract 

QCD sum rules on the light-cone are derived for the sum + /" of the B tt and 
D — > vr form factors taking into account contributions up to twist four. Combining the 
results with the corresponding form factors calculated previously by the same method, 
we obtain the scalar form factors f^. Our sum rule predictions are compared with 
lattice results, current-algebra constraints, and quark-model calculations. Furthermore, 
we calculate decay distributions and the integrated width for the semileptonic decay 
B TTfVr which is sensitive to f^. Finally, the dependence of the sum rules on the 
heavy quark mass and the asymptotic scaling laws are discussed. 
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1 Introduction 



The weak transition B ^ n plays an exceptional role in B physics, in particular at future B 
factories. The amplitude of this transition is given by the hadronic matrix element 

{n{q) I u^,b I Bip + q)) = 2f+{p')q, + + rip"))?, , (1) 

where p + q and q denote the initial and final state four-momenta, respectively, u'jfj.b is the 
relevant weak vector current, and are the two independent form factors. The form factor 
f^{p^) was calculated in |^ using the technique of QCD sum rules on the light-cone. 
Recently, one has also computed the perturbative QCD corrections to /"*" [@, ^. 

In the present paper we complete the calculation of the matrix element (|T]) by deriving the 
corresponding light-cone sum rule for the sum of form factors /"*" + /". This quantity turns 
out to be a pure higher-twist effect. The leading twist-2 contribution vanishes kinematically. 
From /"•" -|- /~ and /"•", one can construct the scalar form factor 

which determines the matrix element of the divergence of the weak vector current: 

(7r(g) I d^{u^,h) I B{p + q)) = (m| - ml)f{p') . (3) 

Using the results on and /° we predict the momentum transfer and lepton energy distri- 
butions as well as the width of the semileptonic decay B nfUr- As a by-product, we also 
obtain the analogous form factors of the D tt transition. 

Furthermore, we investigate the heavy-mass dependence of heavy-to-light form factors. 
The asymptotic scaling laws are determined and found to differ at small and large momentum 
transfer. The origin of this difference is explained in detail. We also study the approach to the 
heavy-quark limit numerically and show that it is reached very slowly. Moreover, the behaviour 
beyound the physical 6-quark mass turns out to be very sensitive to the scale dependence of 
the pion wave functions. 

The paper is organized as follows. In sect. 2, we derive the light-cone sum rule for /+ + /~ 
and compare it with the corresponding sum rule for f'^. The numerical analysis of the new 
sum rule and the resulting prediction of the scalar form factor /° is presented in sect. 3. Sect. 
4 is devoted to the semileptonic decay B nfUr, and sect. 5 to the heavy- mass dependence 
of the form factors. Our conclusions are summarized in sect. 6. 

2 Light-cone sum rule for /+ + f~ 

In order to obtain the QCD sum rule for the form factor combination /+ + appearing in 
(|I]), we follow the method applied to and explained in detail in [JI], ^. The main object of 
investigation is the vacuum-pion correlation function 

F,{p,q) = ^Jd'xe'P^{7c{q)\T{u{x)-f,b{x),m^l5dm\0) 

= F{p\ (p + qf)q^ + F{p\ {p + qf)p^ . (4) 
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Insertion of a complete set of hadronic states with B-meson quantum numbers between the 
currents in (^) entails relations between the physical form factors /"*" and + /~ and the 
invariant amplitudes F and F, respectively. More definitely, for F one finds 



mh{Tn% ~ [p + qy) J s - [p + qy 



"0 



where the term proportional to + / arises from the the contribution of the ground state 
B meson, while the integral over the spectral density represents the contributions from 
excited resonances and continuum states above the threshold energy \fs^. In deriving this 
hadronic representation of F we have used the matrix element (Q) and 

{B I 6i75(i I 0) = m%fB/mb , (6) 

fs being the B meson decay constant. 

In 0, the invariant amplitude F of the same correlation function is calculated 
by expanding the T-product of the currents near the light-cone at x"^ = 0. The leading 
contribution to the operator product expansion (OPE) is obtained by contracting the 6-quark 
fields in (H) and inserting the free 6-quark propagator 

(0|r{W}|0)^./|ie-l±^. (T) 

Substitution of (0) in (^ yields 

FAv.^) = i\ (2^)4(^2 _ p) e-(^-^)-(^m,(vr(g)|..(x)7,75^(0)|0) 

+r(7r(g)|M(x)7^7,75rf(0)|0)^ . (8) 

This approximation is valid in the region of momenta (p + g')^ <^ vn^ and 

p^ <ml- 2mbX , (9) 

X being a mb-independent scale of order Aqcd- Since the pion is on-shell, = vanishes in 
the chiral limit adopted throughout this calculation. The above restrictions ensure that the 
b quark is sufficiently off-shell, and that the resonances in the ub channel are sufficiently far 
away. 

The bilocal vacuum-to-pion matrix elements of light-quark fields encountered on the r.h.s. 
of (|]) are expanded around = leading to a series of contributions with increasing twist. 
The coefficient functions of this expansion can be parametrized by pion wave functions on the 
light-cone |^, 0. Including terms up to order x^, the light-cone expansion of the first matrix 
element in (H) reads 

1 

(7r(g)|n(x)7M75rf(0)|0) = -iq,U J due'^'^^ (<^^(n) + xV(w)) 



+U (x, - ^\ C due'-^^g^iu) . (10) 
V qx I Jo 
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Here, Lp-j^ is the leading twist 2 wave function, while gi and g2 are twist 4 wave functions. 
Upon substitution of 7^7^, = —ia^y + g^^, the second term in (|) is decomposed into the 
matrix elements 

(7r(g) I u{x)i^,d{Q) I 0) = A/i^ du e»''>p(M) (11) 



and 

(7r(g) I M(x)a^,75C^(0) | 0) = i{q,Xy - q.x,)^ £ du e»«>.(M) , (12) 

with /i^ = m1/{mu + md)- In leading order, these matrix elements involve the twist 3 wave 
functions ipp and (pa- It is worth pointing out that the path-ordered gluon operator 



IlG = Pexp{igs dax^A^{ax)} (13) 




ensuring gauge invariance of the above matrix elements is unity in the light-cone gauge, Xf^A^^ = 
0, assumed here. Therefore, the factor Uq is not shown explicitly in (p!0D-(p!2D. 

Substitution of (p!0|)-([T^) in (|^), integration over x and k, and collection of all terms 
proportional to yield the following expression for the invariant amplitude F: 

1 

du \ /i^v9^(n) 



Fqcd{p\ {P + qf) = Uj . _ . \ f^^Mu) + 



^ mf — {p + uqY \ Qu 

xfl ]+ ^"^^^^(^^ |. (14) 

Y ml — {p + uqY j ml — {p + uqyj 

The index QCD has been added to distinguish the above representation of the invariant 
function F in terms of quark and gluon degrees of freedom from the hadronic representation 
given in (^). Note that the twist 2 and 4 wave functions (p-,^ and gi, respectively, do not 
contribute to F . This is obvious from the definition ([T0|). In general, the correlation function 
(^) also receives contributions from gluon emission by the h quark. This correction involves 
quark-antiquark-gluon wave functions as described in 0, . However, direct calculation shows 
that up to the twist 3 and 4 the three-particle correction vanishes in the invariant amplitude F. 
Hence, to twist 4 accuracy, the result for F turns out to be remarkably simple, at least when 
compared with the corresponding expression for the invariant amplitude F given in ||l|, 0] . 

The equality of the two representations (|]) and (0) of F implies a sum rule for + /~ 
which, however, is only useful if one can remove the unknown contributions from the excited 
and continuum states. This is possible to a reasonable approximation by making use of quark- 
hadron duality. Following the standard procedure, the integral in (^ over the hadronic spectral 
function above the ground state is replaced by the corresponding integral over the imaginary 
part of Fqcd- Formally, one can substitute 

~p\p\ s)e{s - s',) = -Im Fqcd{p', s)e{s - s^) , (15) 

TT 

where Sq is an effective threshold parameter separating the duality interval of the ground state 
from the one of the higher states. With this approximation, it is straightforward to subtract 
the contribution of the excited and continuum states from the basic equation given by (Bf) and 
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(|T^. After performing the obligatory Borel transformation in {p + q)"^, one finally arrives at 
the sum rule 



nib 



®0 



7rm 



B 



2 .X / 



lm.FQCD{,P ,s)exp 



M2 



ds 



(16) 



being the Borel mass parameter. 

The remaining task is then to derive ImFQCDip^ ■, s) from (|T^). This is explained below. 
Using {p + uqY = (1 — + u{j) + qY and changing variable from m to s = {ml — p^)/u + p'^ 
one can rewrite (|1^) as follows: 



where 



FQCD{v\{p + q?)= E 



- {p + qfY ' 



Pl(p^s 

2 /t 



S — 



Qu 



ml — p^ 



(mfo - p ) + 2mbg2{u) 



6u 



(17) 



:i8) 



(19) 



The term i = 1 in (|T^ already has the form of a dispersion integral in the variable {p + q)'^- In 
order to achieve this also for the term i = 2 one has to perform a partial integration yielding 
in total: 



FqcDipMp + q)'') 



ds 



s - {p + qY 



Pl(p^s) + 



dp2{p'^,s) 
ds 



°°ds—^ P2(p^s) 



ds \s — {p + qY ) 



(20) 



Since the wave functions ip^j and g2 vanish at u = and u = 1, that is s = cxd and s = ml, 
respectively, as can be seen from the explicit expressions given in the subsequent section, the 
second integral in (pO|) is zero. Hence, the imaginary part of Fqcd can be directly read off 
from the integrand of the first integral: 



It F ^ 2 N ^ 2 N , dp2{p'^,s) 

-ImFQcoiP , s) = pi{p , s) + 

TT ds 



(21) 



Substitution of (|21|) in (|T6D yield; 



rrib 



*0 



m 



B 



dp2{p^,s) 



mh — s 



fBitip') + rip')) = ^ / ( Pi(/, s) + ^^^^^^ ) exp [ 1 ds 



ds 



M2 



ml 



Pl(p^s) + 



P2(p^g) 

M2 



exp 



m. 



M2 



ds + p2{p'^ , Sq) exp 



'_B__fO_ 

M2 



(22) 



In previous applications of QCD light-cone sum rules with the exception of the recent calcula- 
tion of the B ^ p form factors in M, surface terms similar to the last term on the r.h.s. of (E^) 



have been neglected. They originate from higher twist contributions to IirFqcd and play a 
minor role numerically. Nevertheless, in order to subtract the contributions from excited and 
continuum states in the duality approximation consistently we take these terms into account 
in the present calculation. 

The final sum rule for /"*" + /~ follows from (|I8D, (|T9|) and (p2D after returning to the 
variable u: 



/B(r(p') + r(j>')) = %^exp 



m 



B 



+ exp 



6u 




ml — 



u 



.B 



M2 



1 - 



6A 




(23) 



with A = [ml — P^)/{sq — p^). For comparison and later use, we also quote the analogous 
sum rule for /+ obtained in pi 0]: 



/b/+(/ 



2ml 



exp 




ml — j92(l — u) 



+ 



X ( V9^(u) + — 



1 + 



ml + p 
2uM^ 



g2{v)dv 1 



ml + p"^' 



+ t^s^,p',M-^ 



imlgi^u) 



(24) 



Here, we have added the surface term t"*" which was neglected previously, and denoted the 
contribution from the quark- antiquark-gluon wave functions of twist 3 and 4 by Jq. The 
explicit expressions for and Jq can be found in the Appendix. 

Since very recently, the perturbative 0{as) correction to the leading twist 2 piece of the 
light-cone sum rule (|2^) for /"*" is also known |^. However, the corresponding QCD correc- 
tions to the twist 3 term in (p^ ) as well as in the sum rule (^) for + /~ still remain to be 
calculated. Hence, for consistency, we will not include the 0{as) effects in /+ in the present 
analysis. 



3 Numerical results 

For the numerical analysis of the new sum rule (^) we use the same input as in the evaluation 
of the sum rule (p^ in 0]. From experiment we take /,r = 132 MeV and ms = 5.279 GeV, 
whereas the parameters m^ = 4.7 ±0.1 GeV, = 35 =F 2 GeV^ and fs = 140 =F 30 MeV 
are extracted from the QCD sum rule for the correlator of two &75M currents. For consistency, 
the 0{as) correction is not included in the latter two-point sum rule. This is reflected by the 
low value of fs- Because of cancellations of QCD corrections in the ratios of (|^) and fs, 
respectively, (0) and fs, the remaining corrections to the form factors themselves may in fact 
be small. This is precisely what happens in the case of as has been shown in 0, |^. 
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Furthermore, the exphcit expressions for the pion wave functions up to twist 4 are collected 
in 0. Those entering the new sum rule (^) are given below for completeness: 



V9^(M,/i) 
V9p(M,/i) 



6uu 



1 + a2(/i)^[5(M - - 1] + a^{fi)^[21{u - u)^ - U{u - uf + 1]] ,(25) 
Z o J 



l + i?2(y^)2 



3(u-u?-l 



1 r 



+ 34(12)- 35{u - uy - 30{u - uy + 3 , (26) 



6uu 



1 + C^ifi)^ (5(m - uy - 1) + CM^{21{u - uf 



- lA{u-uf + l] 

5 1 6 

-(52(;u)mV + -e(/i)(52(/i)[MM(2 + 13mm) + IOm^ lnM(2 - 3m + -m^) 
2 2 5 



+ 10# lnu(2 - 3m 



6 



-M' 



10 



5^{lJi)uu{u — u) , 



(27) 

(28) 
(29) 



with M = 1 — M and being the renormalization scale. For a detailed discussion of the wave 
functions we refer the reader to the original literature f^, . Recent reviews and references 
can be found in [1^]. Here, the specification of the various coefficients together with a 
few comments may suffice. The terms proportional to the coefficients a^, Bi and Ci represent 
scale-dependent nonasymptotic corrections. They vanish logarithmically as yU ^ 00. 

In leading-order, the value of the scale /i is ambiguous. As a reasonable choice we take 
2.4 GeV. With this choice, the estimate in gives a2{fib) = 0.35, 



mt, - mt 



fJ-b = 

aM = 0.18, B^ifib) = 0.29, B,{fih) = 0.58, C^ifib) = 0.059, CM = 0.034, = 
0.17 GeV^, and e(/ib) = 0.36. Furthermore, the PCAC relation fluT^i^n) = — 2(gg)(/i) and 
{qq){fJ'b) = (—260 ± 10 MeV)'^ can be used to fix the final parameter appearing in (p3[) and 
(^, namely /i^ = fx^ifib) = 2.0 ± 0.25 GeV. 

In the case of the sum rule (p^ for the acceptable range of values of the Borel parameter 
was found to be 8 < < 12 GeV^ In the sum rule (||) for /+ + /" we take the 
same interval after having checked that in this range of the twist 4 contribution does not 
exceed 10%, and that the excited and continuum states (|^) do not contribute more than 30%. 
These are the usual conditions posed in sum rule applications. 

Having specified the necessary numerical input, we are now ready to present quantitative 
results. In Fig. 1, the sum rule ( [23| ) is plotted as a function of Borel parameter M^. One can 
see that the variation is very moderate in the range 8 < < 12 GeV^, at least for < 15 
GeV^. This also applies to the sum rule (^4]). However, at p'^ > 17 GeV^ the dependence 
on becomes strong, indicating that one is getting too close to the physical states in this 
channel. Fig. 2 shows the momentum dependence of the form factors f~^{p'^) + /~(p^), f~^{p'^), 
and /°(p^) for the central value = 10 GeV^. Here, the scalar form factor /° is calculated 
from the other two form factors using the relation (^). In particular, at zero momentum 
transfer we predict 



/+(0) = /0(0) = 0.30, 
/+(0) + r(0) = 0.06. 



(30) 
(31) 



One of the greatest virtues of the sum rule approach is the possibility to estimate the 
theoretical uncertainties in the predictions, at least in principle. In practice, this task is not 
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straightforward and, hence, the estimates should be considered with caution. The uncertainties 
in p^ ) and (|23| ) induced by the input parameters M^, m?,, Sq , fs, and fi^, can be investigated 
by varying these parameters simuhaneously in the sum rules for the form factors and for fs- 
For instance, if the Borel mass is allowed to vary within the interval quoted above, /° is 
found to deviate by ±3% (±5%) at small (large) p"^ from the value obtained with the nominal 
choice = 10 GeV^. The corresponding uncertainty in + /~ can be inferred from Fig. 
1. Furthermore, if and Sq are varied in a correlated way within the ranges given at the 
beginning of this section such that one achieves maximum stability of the sum rule for fs, /° 
changes by about ±2% relative to the value obtained with the nominal choice of parameters. 
This uncertainty is p^-independent. 

Another source of uncertainty is the precise shape of the pion wave functions. Keeping all 
other parameters fixed, we have studied the sensitivity of our results to the nonasymptotic 
terms in the twist 2 and 3 wave functions given in P3|)-(P7D. In Fig. 3 we compare the 
prediction on with the coefficients a^, Bi, and Ci {i = 2,4) as given earlier in this section 
to the result obtained by putting them to zero. The shifts are momentum dependent reaching 
-10 % at = and +5 % at 15 GeV^. The real uncertainty is certainly less than that. 

In addition, there is some uncertainty due to the truncation of the light-cone expansion or, 
in other words, due to the neglect of terms with twist larger than 4. An upper limit on this 
uncertainty should be set by the size of the twist 4 contribution. The influence of the latter 
on the scalar form factor /° is displayed in Fig. 4. The twist 4 terms decrease /° by 2 % at 
small and by 5 % at large p"^. 

Finally, the present lack of knowledge of the perturbative QCD corrections to the twist 3 
contributions in the sum rules (|2^) and (^3]) gives rise to uncertainties which will eventually be 
removed in near future. The scalar form factor /° is concerned in particular at large momentum 
transfer where twist 3 should dominate as expected from the relation (|]). Conversely, at small 
p"^ P coincides with the form factor which receives the leading contribution from the twist 
2 wave function Lp^^. The 0{as) corrections to this piece are known. They change the lowest 
order estimate (^) by only 10 % 

/+(0) = /0(0) = 0.27 . (32) 

As already pointed out, this is due to a remarkable cancellation of the corrections to the sum 
rule (^) for fs and the sum rule for fs in their ratio. Whether or not a similar cancellation 
takes place at the twist 3 level ist questionable and can only be decided by direct calculation. 
For the time being the total uncertainty in /° is estimated to be about 20 % at small p^ and 
30 % at large p^. 

The sum rules (p3|) and (^) for the B —>■ n form factors are easily converted into sum 
rules for the corresponding D —>■ tt form factors. Formally, one only has to replace the flavour 
indices b by c, and B by D. Because of the relatively light charm mass, the region of validity 
of the sum rules covers only the low momentum region < < 1.0 GeV^ of the kinematically 
allowed range of p"^. The values of the input parameters are rric = 1.3 ± 0.1 GeV, = 6 =F 1 

GeV^, and fo = 170 ± 20 MeV. The scale is taken to be fj,c = \/'n^D — ^1 = 1-3 GeV, 
while the fiducial range of the Borel mass is 3 GeV^ < < 5 GeV^. Correspondingly, the 
value of /i^ is lowered to /i7r(/^c) = 1-8 ± 0.5 GeV. The numerical values of the nonasymptotic 
coefficients in the pion wave functions at the scale /Xc are given in 0, There, one can also 
find further remarks on the above choice of parameters. Our numerical predictions for the 
D TT form factors are illustrated in Fig. 5 using the central values for the various input 
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parameters. At zero momentum transfer, we estimate 



/+(0) = /0(0)=0.68, 
r(0) + r(0) = 0.52. 



(33) 
(34) 



4 Predictions for B 



TTTU 



T 



The B ^ 71 form factors can be measured most directly in the weak semileptonic decays 
B — »• 7rli>i where Z = e, /i or r. The distribution of the momentum transfer squared in these 
decays is given by 



dr _ Gl\Vub\H p' 

dp^ 



m. 



fyjm - ml 



247r3 



1 + 



m. 



2p2 



+ 



(35) 



with = {m^g + ml—p^) /2mB being the pion energy in the B rest frame. Another interesting 
observable is the distribution of the charged lepton energy Ei in the B rest frame: 



dV 
dEi 



G\ I Vub I ^ [Pmax ^ 



rmn \ \_ ^ 

2 4 ' 

-4(p2 + 4Ef ) + ^ (SmsEi - 3p2 + Ami) ' ^ 



m 



B 



m 



B 



tip"' 



2m? 



m 



B 



2m\ + p^ — 2m^ — AmBEi + m^ 



f^{p')r{p' 



m. 



ml 



ip' 



m. 



r{p' 



(36) 



with p 



2 

max 



rriBiEi ± 



mf 



+ O(m^). Although the form factors are calculated in the 



chiral limit, otherwise the finite pion mass is taken into account. In the case of light leptons 
I = e, fi the form factor /° or, equivalently, /~ plays a neghgible role because of the smallness 
of the electron and muon masses. Hence, these decay modes can provide information only on 
the form factor In contrast, the decay B nfUr is also sensitive to the form factor /°. 
We therefore concentrate here on the latter case. 

The sum rule results described in the preceding two sections allow to predict the decay 
spectra in the momentum region < < 17 GeV^. In order to include higher momentum 
transfers and to predict integrated widths one has to find another way to calculate the form 
factors up to the kinematical endpoint = (m^ — m^)^ = 26.4 GeV^. In f^, the single-pole 
approximation 

2x fB*gB*Bn 



4.; 



(37) 



2mB*{l — p^ /m\ 

was used. Since the vector B ground state is only about 50 MeV heavier than the pseudoscalar 
B , the B* pole is very near to the endpoint region. Consequently, at maximum p^ the single- 
pole approximation can be expected to be very good. Moreover, the strong B*B'k coupling 
which determines the normalization of the form factor at large p^ can be calculated from the 
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same correlation function from which the sum rule (24) for /+ at low to intermediate 
values of is derived. To this end one employs a double dispersion relation. The method 
and results are described in and reviewed in Extrapolation of the single-pole model to 
smaller matches quite well with the direct estimate from the light-cone sum rule (|2^ ) at 
intermediate momentum transfer = 15 to 20 GeV^. This provides us with a consistent and 
complete theoretical prediction of /''". 

Unfortunately, it is doubtful that a similar procedure can be applied to the scalar form 
factor /°, because the scalar B ground state is expected to be about 500 MeV heavier than the 
pseudoscalar B. Thus, the scalar B pole may be too distant from the kinematical endpoint 
of the B ^ 71 transition for the single-pole approximation to hold. Nearby nonresonant Bn 
states and excited scalar resonances may give comparable contributions. 

Interestingly, there exists a model-independent constraint on the behaviour of the form 
factor at large ~ , i.e., near the kinematical endpoint. The constraint [TT, O 



is derived from a Callan-Treiman type relation obtained by combining current algebra and 
PCAC: 

lim,w, fip') = fB/U. (38) 
In the following we make use of this bound in order to illustrate the sensitivity of the decay 
spectra in B ^ TifVr to the scalar form factor. 

The form factor /° is extrapolated linearly from the value at = 15 GeV^ where the sum 
rules (^) and ( [2^ ) still hold to the value at p^ ~ m| dictated by (P8|). To be conservative 
we take = 150 to 210 MeV in accordance with recent lattice data |[T^ and with QCD sum 
rule estimates (the latter including the perturbative correction, see e.g. This is shown 

in Fig. 6 together with lattice estimates of /°. Obviously, the lattice data favour the lower 
extrapolation. The distributions of p^ and E^- m B ^ tttUt- resulting from the upper and lower 
bounds on /° are plotted in Fig. 7 and Fig. 8, respectively. This study demonstrates that 
measurements of these decay spectra at B factories should provide interesting information on 
the elusive form factor /°. 

For the integrated partial width we predict 

r(5° ^ TT-T+Ur) = 5.7 to 6.5 iKbr Ps'^ , (39) 

where the range corresponds to the two extrapolations considered in Figs. 7 and 8. The 
theoretical uncertainties in the sum rule calculations discussed in sect. 3 are not included. 
The latter drop out to a large extent in the ratio 

^1"^°^"'"^^^ = 0.75 to 0.85 . (40) 

It should be noted that only the numerator is influenced by the scalar form factor. 

The form factor also plays an important role in nonleptonic B decays where it enters 
the factorized two-body amplitudes for B irh. Depending on the mass of the meson h these 
decays probe in the range < p^ < m^. This is similar in D decays. However, there 
the form factor /° cannot be measured independently in semileptonic decays because only the 
electron and muon modes are kinematically accessible. 



5 Dependence on the heavy quark mass 



The light-cone sum rules (^) and (24) offer the possibility to systematically investigate the 
dependence of heavy-to-light form factors on the heavy quark mass. Using the familiar scaling 
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relations for mass parameters and decay constants, to wit 



ms = TTLh + A 



.B 



f. 



B 



= Ib/ , 



= 2mbT 



(41) 
(42) 



where in the heavy quark hmit A, ujq, t, Jb are mft-independent quantities, it is rather straight- 



forward to expand the sum rules in inverse powers of rrib. At p 
given by 



and 



r(o) + r(o) 




0, the leading terms are 



2t 



(43) 



dp exp 



0(m, 



__P_ 

-5/2x 



(44) 



where v'a(l) stands for the derivative dipa/du at the endpoint u = 1. It is important to note 
that whereas the twist 2 and 3 two-particle wave functions survive in the asymptotic limit 
(^31) and (0), the higher-twist and three-particle components are suppressed by one additional 
power of m;,, and die out. 

The asymptotic scaling laws derived above can be understood as follows. At = 0, the 
integration region A = m^/sf < M < 1 in dH) and (H) is rather narrow. In fact, it vanishes 
in the infinite mass limit as 1 — ml/ Sq ~ 2ujQ/mb. In this limit, the asymptotic twist 2 and 
twist 3 wave functions behave like 



iPa ^ {1 — u) Uo/mf, and ipp 



(45) 



Taking into account the extra factor l/m^ multiplying (pp, and noticing that the factors 1/m;, 
times bracket multiplying (p„ in the sum rules approach unity at mi, — >■ oo, one sees that the 
twist 2 and 3 terms lead to the same asymptotic scaling behaviour. The latter is determined 
by a factor m^^ from the integrand, a factor m^^ from the integration region, and a factor m 



1/2 



from 1//b. The fact that the light-cone sum rule predicts f~^{0) 



m, 



'3/2 



was first noticed in 
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The situation changes drastically when the momentum transfer becomes large of order mf. 



To be definite, at the boundary = ml — 2mi,x considered in d), the integration region 
in (23) and (|2^) is finite and independent of m^. Therefore, the asymptotic scaling laws are 
simply determined by the factors in front of the duality integrals in the sum rules, that is 



1// 



B 



m 



1/2 



in the case of and l/m^fB 



V. 



^ for + f~ ■ Explicitly, one obtains 

exp — 



B 



\ f f^du 


X(l - u) 






TU 



X 



6UT 



ifa{u) 



9i{u) - 92{v)dv 
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H — exp 
X 



T 

+ 



1 




92{v)dv 



du 



and 



-1/2 



77 



exp 







\ f /-Mm 








TU 



(46) 



^p(m) + 



6u 



H — exp 



(-J(x + -o)^.(A) 




Ofm, 



-3/2x 



(47) 



where A = + ^^o)- In contrast to the heavy quark hmit at = 0, here also twist 4 

contributes asymptotically. However, the contributions from three-particle wave functions are 
still suppressed by an extra power of l/m^. 

Finally, using the relation (0) it is easy to check that the asymptotic scaling law of the form 
factor coincides with the expectation (|43|) for /"*" at small p^, and with (|47|) for /"*" + /~ at 
large p^. 

The above analysis shows that the light-cone expansion in terms of wave functions with 
increasing twist is consistent with the heavy mass expansion. The higher-twist contributions 
either scale with the same power of rrif, as the leading-twist term, or they are suppressed 
by extra powers of rrib. The sum rules nicely reproduce the asymptotic dependence of the 
form factors on the heavy quark mass as derived in |[Tl|, |1^] for small pion momentum in 
the rest frame of the B meson. In addition, they also allow to investigate the case of large 
pion momentum where neither HQET nor the single-pole model can be trusted. Cleary, as 

— s> excited and continuum states are expected to become more and more important thus 
leading to a break-down of the single-pole approximation. The change in the asymptotic mass 
dependence of the light-cone sum rules when going from large to small momentum transfers 
can be considered as a signal of this break-down. Claims in the literature which differ from 
(^3]) and (^^ are often based on the pole model and therefore incorrect in our opinion. The 
above conclusions corroborate similar analyses carried out for the B K* []TB[ and B ^ p 
P transition form factors. 

In order to clarify the relevance of the asymptotic scaling behaviour in the mass range 



between rric and rub, we have studied the functional dependence of the sum rules (p3D and (p^ ) 
on the heavy quark mass rriQ numerically. The dependence of the parameters rriB, Sq and 
on niQ is described approximately by the relations (|iT| ) using A = 0.6 GeV, ujq = 1.4 GeV, 
r = 1.1 GeV. Together with mq = 4.7 GeV this choice reproduces the central values of etc. 
given at the beginning of sect. 3. In order to consistently include the deviations of the decay 
constant fs from the asymptotic scaling law (^) we have substituted Jb by the corresponding 
two-point sum rule using again the relations ( PT| ) analogously to the procedure described above. 
The logarithmic mass dependence of the wave functions and vacuum condensates through the 
scale hq = ^2mQA is taken into account. Fig. 9 shows the form factor /'*'(0) multiplied 

by the leading power itlq'^ as a function of l/mq. Even at mq > rrib there is still no sign 
that one is approaching the asymptotic limit. On the contrary, the mass dependence of the 
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nonasymptotic terms in the pion wave function becomes more and more important, at least 
for an intermediate mass range. At mg < rrif,, the integration region is still big enough to 
wash out these effects since the wave functions are normalized. In the region between rric and 
rrib, the mass-dependence can be fitted to the following quadratic polynomial in I/ttiq: 

na),nf = 3.3GeVV^ fl - M£!X + , (48) 

indicating the existence of large l/mg corrections in the physical mass range. Similar results 



have been obtained in the case of B ^ K* [16| and B ^ p H form factors 



6 Conclusion 

In this paper, we have derived a new QCD sum rule for the combination /"""(p^) + of 
B ^ TT (and D ^ n) form factors. The light-cone approach used is designed for heavy-to- 
light transitions and incorporates the nonperturbative dynamics in terms of light-cone wave 
functions of the pion. The sum rule is essentially determined by the twist 3 q-q wave functions. 
Terms involving the leading twist 2 wave function and the q-q-g wave function of twist 3 and 
4 are absent, while the twist 4 q-q wave functions give only small contributions. Higher- 
twist components are neglected. The sum rule is valid in the range of momentum transfers 
< p2 < 17 GeV^. 

Combining the new result on + /~ with the corresponding calculation of /"^ in ^ we 
have been able to predict the scalar form factor /°. This prediction is compared with recent 



lattice results [T^ in Fig. 6 and with quark model and different sum rule [|18| estimates 
in Fig. 10. Within the inherent uncertainties of both approaches there is agreement with the 
lattice results. However, the latter tend to be systematically lower than our prediction. This 
could be an indication for the presence of perturbative QCD corrections which still need to 
be calculated. Our result on /° also agrees with the lattice-constrained parametrization of 
this form factor (the pole variant) discussed in The quark model estimate makes use of 
the single-pole approximation /°(p^) = /~^(0)/(l — rn^/p^) with mo = 6.0 GeV. In |12[ it is 



suggested to use instead the relation (^) in order to normalize /° at maximum momentum 
transfer. Extrapolation to intermediate and small then leads to a result very similar to the 
one obtained in fl^ and shown in the figure. Despite of the rough agreement with our sum 
rule result we doubt the validity of the single-pole model at small for reasons explained 
in sect. 5. Finally, in the framework of the heavy-quark-effective-theory one has derived a 
three-point sum rule for /° giving a result in the region < < 10 GeV^ which is about 
30% lower than the expectation from the light-cone sum rule. Concerning similar applications 
of the light-cone sum rules, one should also mention the study of the form factor /~ of the 
B ^ K transition in . 

It would be very interesting to confront these predictions with experimental data, not 
only to test the theoretical methods of calculating form factors, but also since /'^ enters the 
factorized amplitudes for a class of nonleptonic two-body decays. Yet, direct measurements 
of are only feasible in the semileptonic decay B^ — »• 7r~r^z/^. This mode may get in 
experimental reach at future B factories. We have presented the expected decay spectra and 
demonstrated the sensitivity to /°. 

Last but not least, light-cone sum rules for heavy-to- light form factors such as /° provide 
very flexible tools to study the transition from small to large momentum transfers, and to 
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relate the physics of D and B mesons. Moreover, they provide new insights in the heavy 
quark mass dependence of weak matrix elements. 
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Appendix 

Here, we give the explicit expressions for the surface term t+ and the contribution J'q from 
the quark-antiquark-gluon wave functions in the light-cone sum rule 



4ml 



[ml — p'^y 



Ami 



dgi{A) 



M2 



+ 



ml — 



ml — 



M2 



{sq — p'^){ml — p^) du 
2{ml+p^) 



92(A) , (Al) 



JO J («! + was ) 



ml — p^{l — «! — ua^) 
(ai + uaz)M^ 



j 2/3,r / X 



1 - 



2 2 

mi — p 



(A2) 



with "Dui = daida2da36{l — a\ — ai — as). The definitions and functional forms of the twist 
3 wave function 993^ and the twist 4 wave functions (p^ and (p\\ can be found in 0. 
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Figure 1: Form factor (/"^ + / ) (is a function of the Borel parameter at various values of 
the momentum transfer: = (solid), = 10 GeV^ (long- dashed) and p^ = 16 GeV^ 
(short-dashed). 
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Figure 3: Sensitivity of the form factor /° to the light-cone wave functions: nonasymptotic 
corrections included (solid) and purely asymptotic w.f. (dashed) 
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Figure 4: B ^ tt form factor /°; twist 2 and 3 contributions (solid), twist 4 contribution 
(dashed). 
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Figure 6: The B ^ n form factor : direct sum rule estimate (solid) and linear extrapolations 
to the limit ( ^^ ) (dashed). The lattice results are from / |7^/ . 
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Figure 7: Distribution of the momentum transfer squared in B ^ ttti't. The two curves 
correspond to the two extrapolations of /° shown in Fig. 6. 




Figure 8: Distribution of the r-lepton energy in B ^ nfUr- The two curves correspond to the 
two extrapolations of f^ shown in Fig. 6. 
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Figure 9: Dependence of the form factor /^(O) on the heavy quark mass mq for purely 
asymptotic wave functions (solid), and the nonasymptotic corrections included at the scale 

Hq = J2mQA (dashed). 
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